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Abstract. The space-time evolution of J/tp production in central Au+Au collisions 
at RHIC energy is investigated in a transport model. Both gluon dissociation and 
continuous regeneration of J/ips inside the deconfined state are considered. 
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J ft/} suppression [T] is widely accepted as an essential signal of the quark-gluon 
plasma (QGP) formed in relativistic heavy ion collisions. However, the mechanism of 



"^t" ! J 1 4> production in hot and dense medium is still an open question. Different from the 
J/ifj suppression observed at SPS (2] where almost all the charmonia are produced in 
the initial stage via hard processes and then suffer from the anomalous suppression in 
the QGP, there is a remarkable number of charm quarks in the QGP produced at RHIC 
Q\ \ energy and the recombination of those uncorrelated charm quarks offers another origin 

of J ftp production [3]. Both sudden generation on the hadronization surface in grand 
canonical ensemble |4] or canonical ensemble [5] and continuous regeneration 0,[7] inside 
^ ' the QGP are discussed to describe the charmonia production. 

Since the charmonia are so heavy, they are difficult to be fully thermalized in the 
QGP with light quarks and gluons as constitutes, and a natural way to describe the 
charmonium motion in hot and dense medium is through a transport approach. In this 
paper, we investigate the space-time evolution of charmonium production and calculate 
the nuclear modification factor Raa and averaged transverse momentum square (pf) 
for J/ips at RHIC energy in a transport model [8], where the QGP is characterized 
by hydrodynamic equations and the charmonium motion is controlled by a classical 
transport equation. Both initial production and regeneration of charmonia and both 
nuclear absorption and anomalous suppression can be self-consistently considered in the 
model. The leakage effect [H [9j [10], which is especially important for describing the 
transverse momentum saturation at SPS energy [H], is reflected in the free streaming 
term of the transport equation. 

Since a charmonium mass is much larger than the typical temperature of the fireball 
created at RHIC, it is a good approximation to describe the charmonium distribution 
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function fe(p t , x t , r|b) for \P = J /if), if)', Xc in the transverse phase space (p t , x t ) at time 
r and fixed impact parameter b by a Boltzmann equation [8], 

dfa/dr + vy -Vfa = -a^f* + (3q,. (1) 

The second term on the left hand side arises from the free-streaming of \& with transverse 



velocity = Pt/ yPt + m |; which leads to the leakage effect and is important to high 
momentum charmonia [UJ. The first and second terms on the right hand side are 
respectively the suppression and regeneration rates of charmonia. The former is usually 
taken from the gluon dissociation process J/ip + g — > c+c [12]. For J /if) the cross section 
reads 

^M^o ^-f 2 , (2) 
(w/ej/^) 5 

where to is the gluon energy relative to J /if), and ej/^ is the J /if) binding energy Such 
a cross section leads to a flat tail of Raa, and can not explain the plateau structure 
in semi-central collisions and the strong suppression in central collisions observed at 
RHIC [T3| IT3]. Generally, the charmonium binding energy decreases with temperature, 
and the cross section (j2J) is no longer valid above some dissociation temperature Tj, at 
which the binding energy becomes zero. In order to take such an effect into account, we 
modify the loss term a (p t , x t , r|b) [8] by a step function 0(Td — T), 

°* = *k / J0kf^ s ^ < T - 9 ( r * - T ) • < 3 > 

where and are the charmonium and gluon energies, W^(s) is the transition 
probability of the gluon dissociation process as a function of s — (py+p g ) 2 , f g (p g , T, u) 
is the gluon thermal distribution, and T(x 4 , r\b) and u(x. t , r|b) are the local temperature 
and velocity of the hot medium. The gain term (3 can be obtained by the detailed 
balance [8]. 

The local temperature T and fluid velocity u, which govern the thermal gluon 
distribution in a and charm quark distribution in (3 and the suppression and regeneration 
region controlled by the two step functions in ([3]), are determined by hydrodynamic 
equations. We assume that the produced partonic plasma reaches local equilibrium at 
time To = 0.6 fm. After that, the plasma evolves according to the 2+1 dimensional 
Bj or ken's hydrodynamic eauations, 

d T E + V ■ M = -(E + p)/r , 
d T M x + V • (Af,v) = -M x /t - 8 xP , 
d T M y + V • (M y v) = -M y /r - d yP , 

d T R + V ■ (Rv) = -R/r (4) 

with the definitions of E = (e + p)^ 2 — p, M = (e + p)7 2 v and R = 771, where 7 is the 
Lorentz factor, and e, p and v are the energy density, pressure and transverse velocity of 
the QGP. To close the equations, we take the equation of state of ideal gases of partons 
and hadrons with a first order phase transition at T c . The initial condition for the 
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hydrodynamics at RHIC is the same as in Ref [5]. Both the J/ip direct production and 
the feed-down from Xc and ip' are considered in our calculation, and the ratio of them 
in the initial collision is taken as 6 : 3 : 1 |15j . 

In following numerical calculations, we take the dissociation temperature for \c and 
ip' as the critical temperature of the deconfmement phase transition, T£ c = Tf = T c . 
Since the binding energy of J/ i/j in hot and dense medium is estimated to be less than 220 
MeV [HUH?], we take 6jm = 150 MeV and the dissociation temperature 7 1 / = 1.92T C 
to fit the experimental data. From the lattice simulations, the critical temperature is 
taken as T c =165 MeV. Extracted from the experimental data at RHIC energy [T8| [T9], 
the charm quark and charmonium production cross sections at central rapidity region 
in nucleon-nucleon collisions are da ( § N /dy\ y= Q = 120 jib and Buda% N /dy\ y=0 =26A, 4.4 
and 13.2 nb for \l/ = J/ip, 4>' and Xc- 
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Figure 1. The J/t/j nuclear modification factor as a function of time for central 
Au+Au collisions at RHIC energy. 

With the known J ftp distribution fj/^ as a function of time and transverse 
coordinate and momentum, one can easily extract the space-time evolution of 
J ftp production and the final state distributions determined at the hadronization 
hypersurface. The time evolution of the nuclear modification factor Raa for central 
Au+Au collisions at RHIC energy is shown in FigJTJ The evolution starts at the 
thermalization time To- At the beginning, the system is very hot and both the 
suppression and regeneration are significant. The initially produced charmonia, which 
are treated as initial condition /*(pt, x t , T |b) [8] of the transport equation ([!]), suffer 
from the strong gluon dissociation and the number drops fast. Since the regeneration 
is proportional to the square of charm quark number but the suppression is linear in 
gluon number, the number of the survived regenerated J/ips increases with time fast 
in the beginning. With the expansion of the system, the temperature decreases, and 
the suppression and regeneration become smooth and finally get saturated at the phase 
transition temperature T c which corresponds to the time r ~ 4 fm/c. Since we did not 
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consider the suppression and regeneration in hadronic phase, the initial production and 
regeneration are time independent after the phase transition. Due to the competition 
between the initial production and regeneration, the total J/ip production is almost a 
constant in the whole evolution. 

Where are the finally observed J/ifjs from, from the central part of the fireball or 
the surrounded region? We calculated the source distribution of the observed J/ips 
from the function fji^ in the limit of r — > oo. The J/ip number density as a function 
of the fireball radius in the transverse plane is shown in FigfSl The initially produced 
charmonia in the central region are fully eaten up by the extremely hot QGP in the initial 
stage when the temperature is larger than T^. Since regeneration happens continuously 
in the whole QGP region, those regenerated J/ifjs in the later stage with temperature 
T <Td have the probability to survive, and therefore the finally observed J/i/js coming 
from the center of the fireball are mostly regenerated. In the periphery of the fireball 
where the temperature is low, both the suppression and regeneration are weak, and the 
observed J/ips from this region are dominated by the initial production. For those J/ips 
coming from the middle part of the fireball, both the initially produced and regenerated 
J ftps are important. 
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Figure 2. The observed J/tp number density as a function of the fireball radius in 
the transverse plane for central Au+Au collisions at RHIC energy. 

We now turn to discuss the centrality dependence of the nuclear modification factor 
Raa and the averaged transverse momentum square (pf) for J/tfis. The calculated Raa 
as a function of the participant nucleon number N p for Au+Au collisions at RHIC energy 
and the comparison with the experimental data [T3J [H] are shown in the left panel of 
Figj3j Since both the suppression and regeneration increase with centrality, the yield of 
the initially produced J/ips drops monotonously and the yield of the regenerated J ftps 
goes up monotonously. For peripheral and semi-central collisions where the temperature 
of the system is low, the J/ip production is controlled by the initial production, while 
for central collisions where the charm quark density is high, both the initial production 
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and regeneration are important. From the data, there exists a flat structure in between 
N p = 60 and 170, and then the Raa decreases further for N p > 170. The plateau can be 
explained by the competition between the initial production and regeneration, and the 
further decrease in large N p region is due to the fact that the maximum temperature 
of the fireball created in central collisions is higher than the dissociation temperature 
T(i- For extremely central collisions around N p — 350, the difference between the data 
and the theory is probably due to the geometry fluctuations which are proved to be 
important at SPS energy [ITJ [20] . 
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Figure 3. The nuclear modification factor Raa (left panel) and averaged transverse 
momentum square (pf) (right panel) as functions of participant nuclcon number N p in 
mid-rapidity region for Au+Au collisions at RHIC energy. The data are from PHENIX 
collaboration [131 114j . 

Different from the global yield, the momentum distribution is more sensitive to the 
mechanism of J/ip production and suppression. Our theoretical calculation of J/ip (pf) 
as a function of centrality for Au+Au collisions at RHIC energy and the comparison with 
the experimental data [13] are given in the right panel of Figj3j Similar to the calculation 
at SPS energy [TT], the (pf) of the initially produced J ftps increases with centrality and 
gets saturated in central collisions, due to the gluon multi-scattering[2T] with nucleons in 
the initial state and the leakage effect which is especially important for high momentum 
J/ips. Since we have assumed kinetic equilibration of charm quarks in the QGP, the 
transverse momentum distribution of the regenerated J / ips is controlled by the charm 
quark thermal distribution, and the corresponding (pf) is only about half of that for 
initially produced ones. The total (pf) is dominated by the competition between the two 
production mechanisms. At high N p the contribution from the regeneration becomes 
important, and the competition leads to a transverse momentum suppression in central 
collisions which agrees well with the data. 

We calculated also the Raa and (pf) for Cu-Cu collisions at RHIC energy. Since 
the colliding energy for Au+Au and Cu+Cu collisions are the same, and the nuclear 
geometry is mainly reflected in the participant nucleon number, the Raa and (pf) for the 
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two kinds of nuclear collisions are almost the same in the centrality region of N p < 110 
where 110 is the maximum N p for Cu+Cu collisions. Our calculation agrees with the 
experimental data [22]. 

In summary, we have calculated the space-time evolution and the observed nuclear 
modification factor and averaged transverse momentum square of J ftps produced in 
relativistic heavy ion collisions at RHIC energy. While the initial production dominates 
the peripheral and semi-central collisions, the initial production and regeneration are 
almost equally important for central collisions. 
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